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"nVhen  brittle  metals  are  subjected  to  blast  loads,  they  tend  to  shatter  Into  a large  number  of  very 
small  particles,  the  number  and  size  of  the  particles  generated  by  such  an  explosive  process  Is  well- 
characterized  by  the  semi -empirical  cumulative  distribution  function  of  Mott.  This  Auction,  however, 
contains  no  Information  regarding  either  the  manner  of  break-up  or  the  reason  for  a particular  distri- 
bution. The  present  experimental  investigation  was  undertaken  to  clarity  the  role  played  by  dynamic 
material  properties  in  the  ultimate  distribution  of  particle  sizes  when  brittle  metals  fracture  under  - — 
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explosive  loads.  The  phenomenology  of  brittle  fragmentation  was  examined  previously,  and  a simple 
model  was  developed  utilizing  the  concept  of  crack  branching.  Using  statistical  arguments  to  establish 
relationships  between  the  energies  associated  with  crack  branching  and  the  Mott  parameters  that 
characterize  the  particle  distribution,  that  model  predicts  that  the  average  mass  of  a particle  from  the 
shattered  body  is  proportional  to  the  fourth  power  of  the  crack  branching  stress  intensity  factor.  To 
examine  the  physical  relationships  between  material  properties  and  particle  size  distributions  from 
blast-loaded  structures,  two  hypereutectoid  steels  (FS-01  and  HF-1)  were  heat-treated  to  a number 
of  brittle  conditions  and  evaluated  in  terms  of  static  tensile  properties,  plane  strain  fracture  toughness, 
crack  branching  behavior,  and  explosive  fragmentation.  The  distribution  parameters  did  not  correlate 
with  strength,  toughness,  or  hardness,  but  the  coarseness  of  the  particles  appeared  to  increase  with 
increasing  tensile  elongation^The  ratio  (K^/Kg;)3-4  was  experimentally  determined  to  be  linearly 
related  to  the  average  p articled  ass  in  reasonably  good  agreement  with  the  model  prediction,  where 
Kh  is  the  dynamic  stress  intensity  factor  for  incipient  microbranching  (hackle). 
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CORRELATIONS  OF  CRACK  BRANCHING  PROPERTIES 
WITH  MOTT  FRAGMENTATION  PARAMETERS 
IN  HYPEREUTECTOID  STEELS 


INTRODUCTION 


For  the  past  several  decades,  the  major  concern  of  analytical  and 
experimental  fracture  mechanics  has  been  the  definition  of  those  condi- 
tions that  precipitate  crack  instability  in  structures.  Much  recent 
interest  has  focused  upon  the  conditions  that  govern  the  propagation 
behavior  of  dynamic  cracks.  The  impetus  for  this  shift  has  been  pro- 
vided by  the  need  to  improve  crack  arrest  properties  of  materials, 
particularly  those  used  in  the  power  and  transportation  industries. 
Another  aspect  of  material  behavior  dominated  by  crack  propagation 
characteristics  rather  than  by  crack  initiation  conditions,  however,  is 
the  susceptibility  of  high-strength  materials  to  shatter  rather  than 
simply  to  yield  or  to  rupture.  Such  phenomena  are  characterized  by 
complex  interactive  networks  of  simultaneously  propagating  fractures 
marked  by  frequent  crack  bifurcation  events. 

N.F.  Mott  considered  the  problem  of  natural  fragmentation  of 
metals  due  to  blast  loading  and  subsequently  proposed  a dynamic  model 
of  the  fracture  of  an  explosively-loaded  ductile  metal  ring  assuming  a 
time-dependent  probability  of  random  fracture  after  a critical  strain 
had  been  reached  in  the  ring.(l)  From  this,  the  form  of  the  particle 
size  distribution  can  be  determined  empirically.  When  the  resultant 
distribution  of  particle  sizes  is  analyzed,  the  frequency  of  occurrence 
of  particular  masses  follows  a cumulative  distribution  described  by  a 
simple  exponential  relationship: 

N “ N exp  (-/m'/vO  (1) 

m o 

where  Nm  - total  number  of  particles  having  mass  greater  than  m,  and 
N0  and  y are  experimentally  determined  parameters  representing  the 
total  number  of  particles  and  half  the  average  mass  of  the  particles. 

This  relationship  describes  the  explosive  break-up  of  ductile 
metals  quite  well.  However,  Mott's  model  is  based  upon  instantaneous 
rupture  of  thin  rings  and  uses  the  root-mean-square  scatter  of  the 
static  fracture  strain  as  the  criterion  for  generating  a particular 
particle  size  distribution.  It  is,  therefore,  somewhat  surprising  to 
find  that  equation  (1}  also  describes  the  distribution  of  particle 

Note:  Manuscript  submitted  May  16,  1979. 
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sizes  that  results  from  the  explosive  loading  of  extremely  brittle 
thick-walled  cylinders.  It  is  a well-known  fact  that  cracks  in  highly 
stressed  brittle  materials  frequently  bifurcate;  therefore,  the  stress 
intensity  factor  for  crack  branching  in  a brittle  material  should  be 
expected  to  be  functionally  related  to  the  statistical  parameters  that 
describe  the  resultant  explosive  break-up  of  the  material. 

On  the  basis  of  preliminary  analyses  of  the  possible  relation- 
ships between  crack  branching  stress  intensity  factor  (Kjg)  and  the 
statistical  Mott  parameters,  Weimer  and  Rogers  (2)  suggested  that  the 
average  particle  mass  after  break-up  should  be  proportional  to  the 
fourth  power  of  Kjg.  They  reported  experimental  results  for  two  hyper- 
eutectoid  steels  which  indicated  the  exponent  to  be  4.55  (+0.17). 
However,  subsequent  experiments  and  corrections  to  the  stress  analysis 
reduced  this  exponent  to  3.03  (+0.65)  with  even  greater  intuitive  un- 
certainty than  this  standard  deviation  would  indicate,  because  fully 
developed  branch  cracks  seldom,  if  ever,  occurred  in  some  of  the 
material  conditions  tested. 

The  objective  of  this  investigation  is  to  clarify  the  role  played 
by  dynamic  fracture  properties  of  a brittle  metal  in  the  distribution 
of  particle  sizes  that  result  when  this  metal  is  subjected  to  blast 
loads.  To  accomplish  this  end,  a simple  model  of  the  explosive  shatter 
of  a metal  cylinder  has  been  constructed  and  described  in  detail  else- 
where. Probabilistic  arguments  were  used  to  demonstrate  that  the 
well-established  Mott  cumulative  distribution  function  could  be  formally 
extracted  from  energy  considerations  of  crack  branching  phenomena.  It 
was  concluded  that  the  average  mass  and  total  number  of  particles  from 
a shattered  body  did  indeed  bear  fourth-power  relationships  to  the  ratio 
of  the  crack  branching  stress  intensity  factor  to  the  plane  strain 
fracture  toughness  of  a given  material.  However,  the  notion  of  what 
constitutes  crack  branching  in  the  explosive  loading  environment  was 
revised  to  include  hackle,  or  incipient  bifurcation,  and  experiments 
were  conducted  to  evaluate  this  approach. 


EXPERIMENTAL  PROCEDURE 


Selection  of  Materials:  Two  hypereutectoid  steels  were  selected 
for  this  investigation.  One  is  FS-01  steel  (0.92  C,  0.30  Si,  1.15  >fci, 
0.50  W,  0.50  Cr,  0.20  V),  a fine-grained  oil  hardening  tool  steel.  The 
other  is  HF-1  steel  (1.08  C,  0.90  Si,  1.88  Mn) , a relatively  large- 
grained silico-manganese  steel.  FS-01  steel  was  selected  because  rather 
simple  thermal  treatments  could  be  expected  to  produce  the  high  strength 
and  hardness  levels  at  which  it  was  anticipated  that  crack  branching 
could  be  easily  and  consistently  produced.  In  addition,  the  nominal 
composition  was  similar  to  that  of  the  steel  used  by  Anthony  (4)  in  his 
crack  branching  investigations.  Because  of  the  paucity  of  experimental 
data  on  crack  branching  in  metals,  it  was  considered  important  at  the 
outset  to  have  some  basis  for  comparison  of  metallurgical  effects. 
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HF-1  steel  was  selected  because  of  extensive  experience  with  particle 
size  distributions  arising  from  explosive  loading.  It  was  also  valu- 
able to  have  established  heat  treatment  schedules  for  producing 
desired  distributions  in  order  to  economize  on  experimental  design. 

Materials  Characterization : The  materials  used  in  this  program 
were  acquired  from  commercial  sources.  The  FS-01  steel  was  available 
both  in  4.8  mm  thick  surface-ground  plate  stock  and  63.5  mm  diameter 
round  bar  stock,  while  the  HF-1  steel  was  available  only  in  4 x 4 RCS 
(10  cm  x 10  cm  round-corner-square)  bar  stock.  The  heat  treatments 
selected  for  this  investigation  were  such  that  the  hardness  values 
obtained  for  FS-01  steel  spanned  the  range  obtained  by  Anthony  and  in 
addition  included  one  condition  of  considerably  greater  ductility  than 
he  utilized.  It  was  also  required  that  the  HF-1  steel  be  heat-treated 
to  conditions  identical  to  those  obtained  in  prior  work  (5)  so  that  the 
explosive  tests  need  not  be  repeated  for  that  material.  The  actual 
selection  was  predicated  on  the  availability  of  fragments  recovered 
from  the  earlier  HF-1  tests  as  well  as  by  the  large  spread  in  size 
distribution  parameters.  The  designations  and  schedules  for  final 
heat  treatments  are  listed  in  Table  1.  These  designations  are  retained 
throughout  this  work  and  are  used  to  identify  material  and  its  condi- 
tion. The  four  FS-01  heat  treatments  produce  very  fine-grained  tem- 
pered martensites  at  several  strength  levels.  The  isothermal  treat- 
ments of  HF-1  cause  proeutectoid  carbides  to  precipitate  at  the  grain 
boundaries.  Moreover,  the  austenitizing  temperatures  are  sufficiently 
high  for  A and  C that  significant  grain  growth  occurs.  The  resultant 
microstructures  are  tempered  martensite  matrices  with  carbide  networks 
along  prior  austenite  grain  boundaries. 

Conventional  engineering  mechanical  properties  were  measured  in 
accordance  with  ASTM  specifications,  including  tensile  properties, ^ 
hardness,  (7)  an(j  fracture  toughness. (8)  The  0.2%  offset  yield  strength, 
ultimate  tensile  strength,  and  percent  elongation  in  a 51  mm  gage 
length  were  measured  using  pin-loaded  flat  tensile  specimens.  Tests 
were  conducted  at  room  temperature  ( -22C)  at  strain  rates  of  10-^  sec-^. 
Hardness  was  measured  using  the  C-scale  of  a standard  Rockwell  Hardness 
Tester.  The  plane  strain  fracture  toughness  was  measured  at  room  tem- 
perature using  single-edge-notch  tension  specimens  that  were  pin-loaded 
quasistatically.  A summary  of  all  mechanical  properties  measured  is 
presented  in  Table  2. 

Explosive  Testing:  Selection  of  explosives  is  based  upon  the 
nature  of  the  task  to  be  performed.  If  the  requirement  is  heat,  de- 
formation, or  simply  lifting  or  heaving  ability,  the  explosive  must 
have  a large  total  work  capacity,  which  is  simply  a high  caloric  out- 
put per  gram  of  explosive.  On  the  other  hand,  when  the  requirement  is 
for  high  shattering  effects,  the  explosive  must  have  a high  rate  of 
delivery  of  that  caloric  output.  The  shattering  power  of  an  explosive 
is  termed  brisance  and  is  linearly  related  to  the  detonation  rate  (9) . 
Since  the  expressed  intent  of  this  work  is  to  investigate  shattering 
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TABLE  1. 


HEAT  TREATMENT  SCHEDULES  AND  DESIGNATIONS 


FS-Ol(A) 


HF-l(A) 


Preheat  to  538C 
Transfer  to  800C,  1 hour 
Marquench 

Temper  at  191C,  2 hours 
Air  Cool 


927C,  1 hour 
Transfer  to  760C, 
Oil  quench 
Temper  at  482C,  1 
Air  Cool 


FS-Ol(B) 


HF-l(B) 


Preheat  to  538C 
Transfer  to  800C,  1 hour 
Marquench 

Temper  at  274C,  2 hours 
Air  Cool 


774C,  1 hour 
Oil  quench 
Temper  at  593C,  1 
Air  Cool 


FS-Ol(C) 


HF-l(C) 


Preheat  to  538C 
Transfer  to  800C,  1 hour 
Marquench 

Temper  at  333C,  2 hours 
Air  Cool 


982C,  1 hour 
Transfer  to  760C, 
Oil  quench 
Temper  at  649C,  1 
Air  Cool 


FS-Ol(D) 


Preheat  to  538C 
Transfer  to  800C,  1 hour 
Marquench 

Temper  at  438C,  2 hours 
Air  Cool 


1 hour 
hour 


hour 


1 hour 
hour 
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(a)  Based  on  at  least  three  tests  and  stated  as  the  mean  value  (+  one  standard  deviation) . 


effects  in  brittle  metals.  Composition  B was  selected  as  the  test  ex- 
plosive. Composition  B is  a derived  or  secondary  high  explosive  con- 
sisting of  64%  cyclotrimethylene  trinitramine  (RDX)  and  36%  trinitro- 
toluene (TNT).  It  has  an  initial  density  of  1.71  g/cm3,  a detonation 
velocity  of  8.01  mm/ps  which  is  supersonic  with  respect  to  the  metal, 
and  a peak  explosion  pressure  of  approximately  30  GPa  (more  than 
4,000,000  psi). 

Three  test  cylinders  were  fabricated  for  each  material  condition 
listed  in  Table  1.  They  were  machined  directly  from  the  bar  stock  to 
dimensions  slightly  greater  than  shown  in  Figure  1 which  shows  the 
geometry  of  a test  cylinder  when  charged  with  explosive.  After  heat 
treatment,  each  cylinder  was  surface  ground  to  final  dimensions  and 
shipped  to  Denver  Research  Institute  where  it  was  loaded  with  explosive, 
detonated,  and  partially  recovered  under  service  contract  DAAA25-72-R- 
0236.  The  C/M  ratio  (charge  mass  to  cylinder  mass)  was  held  constant 
at  approximately  0.3. 

Composition  B is  a cast  explosive.  It  was  extended  beyond  the 
cylinder  ends  to  approximate  more  closely  a plane  wave  when  the  deton- 
ation front  reaches  and  departs  the  cylinder.  The  taper  simply  facili- 
tates removal  of  the  end  cap  used  in  casting. 

A soft  partial  recovery  system  was  used  to  collect  particles  of 
the  shattered  cylinder  as  shown  in  Figure  2.  This  technique  has  evolved 
from  experience  and  direct  comparison  with  total  recovery  systems;  the 
results  are  equivalent  but  simpler  to  obtain.  The  system  is  called  soft 
in  relation  to  recovery  in  sand  or  sawdust.  Basically  it  consists  of  a 
fragment  shield  and  recovery  tank  located  far  enough  from  the  exploding 
cylinder  that  the  charge  can  be  considered  as  a point  source.  In  the 
shield  there  is  a window  that  subtends  an  angle  of  30°  at  the  source 
and  allows  fragments  to  pass  through  to  the  recovery  tank.  The  latter 
is  filled  with  water  and  covered  with  a high  density  urethane  foam  pad 
as  shown. 

Trisodium  phosphate  (NajPO^)  and  sodium  nitrite  (NaN02)  were  added 
to  water  at  a rate  of  1.9  grams  of  each  per  liter  of  water  to  inhibit 
rusting  of  the  recovered  particles.  The  actual  recovery  is  on  a framed 
50  mesh  brass  screen  submerged  in  the  tank  and  covering  the  entire 
bottom. 

The  particles  are  removed  from  the  recovery  screen  by  an  electro- 
magnet, dried,  and  sorted  into  size  groups  by  sieving;  after  counting 
they  are  batch-weighed  and  returned  for  analysis  of  the  distribution 
data.  The  particle  sizes  were  grouped  according  to  standard  U.S.  Series 
sieve  sizes  as  follows: 
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FIGURE  1.  SPECIMEN  GEOMETRY  FOR  EXPLOSIVE  STRAIN  RATE  TESTS 


FIGURE  2.  TEST  SETUP  FOR  SOFT  RECOVERY  OF  MATERIAL 
FROM  EXPLOSIVE  STRAIN  RATE  TESTS 
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Mesh 

Number 


Opening 

(mm) 


3 

6.68 

4 

4.70 

6 

3.33 

8 

2.36 

10 

1.98 

12 

1.65 

16 

1,17 

The  particle  size  distributions  obtained  in  triplicate  tests  on  the 
seven  material  conditions  are  presented  In  Tables  3 and  4.  These  data 
have  been  corrected  for  efficiency  of  recovery  and  normalized  to  full- 
cylinder  recovery,  such  that  the  average  mass  corresponding  to  any  in- 
terval (sieve  size)  is  identical  to  that  calculated  from  the  initial 
partial  recovery  data.  The  actual  parameters  that  characterize  the 
distribution  are  obtained  from  analysis  of  Mott  plots  of  these  data. 

To  construct  a Mott  plot,  sieves  are  ordered  from  coarsest  to 
finest  mesh  and  the  average  mass  of  particles  retained  on  each  sieve  is 
calculated.  The  logarithm  of  the  cumulative  number  of  fragments 
collected,  is  plotted  as  the  ordinate  and  Sm  is  plotted  as  the  abscissa. 
Because  this  distribution  is  highly  skewed,  any  undue  bias  in  assigning 
numbers  to  the  data  points  produces  highly  misleading  results.  Such  a 
bias  results  if  we  take  Nm  to  be  the  total  number  of  particles  accumu- 
lated on  a certain  number  of  sieves  because  it  positions  a data  point 
at  the  edge  of  an  interval  in  a situation  where  the  number  of  counting 
bins  is  quite  limited.  To  correct  this  the  number  of  particles  assigned 
to  the  first  sieve  is  taken  to  be  half  the  number  counted.  This  assumes 
the  distribution  is  linear  over  small  intervals,  and  half  the  particles 
will  have  mass  less  than  the  computed  average  mass  for  that  sieve.  Then 
half  the  particles  on  the  next  smaller  sieve  are  added  to  all  the  parti- 
cles on  the  larger  sieve.  This  number  is  plotted  as  the  average  mass  of 
the  second  sieve.  This  process  is  continued  for  all  sieves  and  yields 
graphs  such  as  presented  in  Figure  3,  which  also  shows  the  particles 
recovered  from  a 30°  sector  of  a typical  shattered  cylinder  in  relation 
to  the  size  of  the  original  cylinder.  The  lines  through  the  data  in  the 
figures  are  the  best  straight  lines  fitted  by  linear  regression  to  the 
function  Nn,  ■ N0e"^o/P , the  Mott  distribution.  Correlations  to  this 
function  averaged  96%  or  better  over  the  seven  material  conditions 
tested.  The  Mott  number  parameter,  N0,  is  simply  the  extrapolated 
intercept  of  these  curves  with  the  ordinate.  The  Mott  size  parameter, 
y,  is  obtained  as  the  square  of  the  reciprocal  slope  of  these  curves. 

A summary  of  Mott  parameters  for  all  tests  is  presented  in  Table  3. 
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TABLE  3 


FS-01  STEEL  PARTICLE  SIZE  DISTRIBUTIONS 


u.s. 

Sieve 

Size 

Particle 

Count 

Ni 

Total 

Mass 

Mi  (g) 

Particle 

Count 

Ni 

Total 

Mass 

Mi  (g) 

Particle 

Count 

Ni 

Total 

Mass 

Mi  (g) 

Cylinder 

A- 1 

Cylinder 

A- 2 

Cylinder  A-3 

3 

0 

0 

0 

0 

0 

0 

4 

108 

62.8 

89 

64.6 

129 

76.8 

6 

1080 

248.7 

800 

210.1 

1010 

243.0 

8 

1884 

194.2 

2074 

208.6 

1903 

197.0 

10 

1260 

60.1 

1644 

73.6 

1722 

83.7 

12 

1488 

42.6 

1852 

50.1 

1087 

29.7 

16 

3540 

49.4 

4030 

55.0 

2485 

33.1 

Cylinder 

B-l 

Cylinder 

B-2 

Cylinder  B-3 

3 

11 

15.5 

29 

30.6 

22 

27.3 

4 

227 

129.3 

218 

98.1 

197 

114.4 

6 

941 

229.7 

988 

239.1 

1064 

265.1 

8 

1818 

185.3 

1874 

197.2 

1470 

150.2 

10 

1082 

50.9 

1191 

53.7 

1272 

55.7 

12 

952 

26.7 

886 

25.8 

899 

24.7 

16 

2337 

32.5 

2208 

28.3 

2161 

28.5 

Cylinder 

C-l 

Cylinder 

C-2 

Cylinder 

C-3 

3 

0 

0.0 

26 

38.3 

29 

38.3 

4 

236 

156.8 

261 

152.5 

315 

208.0 

6 

1095 

271.2 

899 

246.9 

888 

221.5 

8 

1519 

159.5 

1433 

145.3 

1232 

139.8 

10 

930 

45.6 

808 

39.1 

687 

33.0 

12 

577 

15.8 

1003 

27.7 

831 

15.9 

16 

1743 

23.2 

2046 

27.2 

1389 

21.5 

Cylinder  D-l 

Cylinder 

D-2 

Cylinder  D-3 

3 

0 

0 

41 

92.4 

42 

84.9 

4 

474 

407.3 

409 

312.2 

465 

322.6 

6 

536 

153.9 

600 

162.1 

691 

173.0 

8 

686 

71.5 

655 

81.9 

663 

65.1 

10 

337 

17.7 

355 

16.9 

465 

22.0 

12 

362 

11.5 

395 

11.5 

353 

10.2 

16 

873 

25.0 

777 

10.9 

747 

12.0 
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TABLE  4 


HF-1  STEEL  PARTICLE  SIZE  DISTRIBUTIONS 


u.s. 

Sieve 

Size 

Particle 

Count 

Ni 

Total 

Mass 

Mt  (g) 

Particle 

Count 

Ni 

Total 

Mass 

Mi  (g) 

Particle 

Count 

Ni 

Total 
Mass 
Mi  (g: 

Cylinder  A-l 

Cylinder  A-2 

Cylinder  A-3 

3 

0 

0 

0 

0 

0 

0 

4 

136 

76.3 

157 

67.7 

179 

91.2 

6 

829 

232.3 

964 

228.2 

930 

225.6 

8 

1902 

194.6 

2024 

212.7 

2207 

199.7 

10 

1399 

70.9 

1325 

61.0 

1551 

69.6 

12 

1603 

44.5 

1518 

42.0 

1742 

46.6 

16 

3818 

50.1 

3783 

51.8 

3555 

47.9 

Cylinder  B-l 

Cylinder  B-2 

Cylinder  B-3 

3 

50 

85.9 

54 

64 . 6 

63 

79.6 

4 

315 

282.0 

326 

198.9 

403 

260.8 

6 

693 

188.5 

957 

269.3 

755 

206.2 

8 

920 

97.5 

1120 

120.1 

1069 

103.2 

10 

391 

20.0 

554 

26.1 

654 

28.0 

12 

454 

14.2 

511 

15.0 

428 

10.6 

16 

1009 

14.2 

674 

9.8 

755 

9.9 

Cylinder  C-l 

Cylindei 

r C-2 

Cylinder  C-3 

3 

0 

0 

0 

0 

0 

0 

4 

341 

237.8 

389 

250.8 

236 

152.5 

6 

921 

241.6 

663 

175.9 

957 

283.7 

8 

1401 

136.9 

1554 

181.8 

1586 

147.4 

10 

681 

30.5 

834 

40.4 

. 1022 

50.4 

12 

770 

21.8 

663 

18.7 

957 

27.3 

16 

1842 

24.0 

1589 

21.2 

1795 

25.4 

11 


FIGURE  3. 


Mott  Plots  for  FS-01  (A)  Steel.  ^ Is  the  total  number 
of  particles  having  mass  greater  than  m.  Partial  recovery 
data  was  normalized  to  full  recovery.  The  particles 
recovered  from  a 30°  sector  of  a typical  cylinder  are 
shown  to  provide  a size  perspective. 
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TABLE  5 


MOTT  PARAMETERS  FOR  EXPLOSIVE  TESTS 


Heat 

Test 

Initial 

Efficiency 

Mott  Number 

Mott  Size 

Treatment 

Cylinder 

Mass 

of 

Parameter 

Parameter 

Designation 

Number 

Mo  (8) 

Recovery 

No 

W(mg) 

FS-01  (A) 

A-l 

692 

1.000 

20265 

17.4 

A- 2 

701 

0.810 

18011 

20.1 

A- 3 

691 

0.927 

18453 

19.4 

FS-01  (B) 

B-l 

692 

1.109 

15256 

23.1 

B-2 

688 

0.826 

14647 

23.8 

B-3 

688 

1.094 

13562 

25.5 

FS-01  (C) 

C-l 

690 

1.018 

10623 

34.1 

C-2 

695 

0.920 

10615 

33.2 

C-3 

695 

0.837 

9610 

34.9 

FS-01  (D) 

D-l 

695 

0.962 

3999 

104.9 

D-2 

695 

0.880 

4002 

98.3 

D-3 

696 

0.939 

4474 

88.9 

HF-1  (A) 

A-l 

707 

0.736 

14958 

25.3 

A- 2 

701 

0.830 

16528 

22.8 

A- 3 

708 

0.838 

17578 

21.8 

HF-1  (B) 

B-l 

710 

0.793 

5460 

65.5 

B-2 

709 

0.920 

5995 

75.8 

B-3 

707 

0.795 

5771 

70.9 

HF-1  (C) 

C-l 

708 

0.792 

8290 

46.7 

C-2 

706 

0.875 

8565 

44.8 

C-3 

705 

0.763 

9953 

35.0 

13 


Crack  Branching  Studies:  After  completing  the  characterization  of 
the  materials  in  terms  of  both  mechanical  properties  and  particle  size 
distributions,  the  crack  branching  stress  intensity  factors  for  each 
material  were  determined.  Although  it  is  rather  easy  to  induce  crack 
branching  in  a brittle  material,  the  analysis  used  to  calculate  the 
stress  intensity  factor  at  branching  in  all  but  the  simplest  geometries 
is  often  suspect.  The  dynamics  of  fracture  in  finite  specimens  is  not 
well  described;  most  of  the  analyses  are  for  infinite  specimens  or 
physically  artificial  cracks.  One  is  forced  to  construct  an  operational 
definition  of  the  crack  branching  stress  intensity  factor  for  finite 
specimen  geometry.  The  most  practical  approach  to  inducing  branch 
cracks  in  a metal  is  to  load  a blunt  notched  specimen  in  tension.  Un- 
fortunately, acceleration  of  the  resulting  crack  is  far  more  intense 
than  would  be  experienced  by  a sharp  natural  crack,  and  the  average 
strain  energy  density  in  the  material  is  considerably  higher.  Changing 
the  notch  stress  concentration  factor  changes  the  fracture  load  but  pre- 
sumably not  the  crack  branching  stress  intensity  factor. 

The  specimen  developed  for  these  studies  is  a modified  single-edge- 
notch  (SEN)  tension  specimen  with  a blunt  notch  of  constant  root  radius 
equal  to  1.27  mm.  This  specimen  and  the  formulation  of  the  static  stress 
intensity  factor  for  arbitrary  crack  length  are  illustrated  in  Figure  4. 

The  high  hardness  of  the  materials  being  investigated  creates  a need 
for  special  gripping  arrangements  to  load  a plate  in  tension.  Anthony 
(4)  solved  this  by  bonding  soft  mild  steel  plates  to  the  specimen  grip 
areas  with  an  epoxy  resin.  The  approach  in  the  current  work  has  been  to 
use  pin-loading,  reducing  the  gage  section  sufficiently  for  large  frac- 
ture loads  to  be  applied  without  fracturing  the  specimen  at  the  pin  holes. 
The  length  of  the  reduced  section,  i.e.,  the  gage  length,  was  20  cm  and 
the  separation  of  load  points  was  30  cm.  This  assures  that  the  tensile 
field  in  the  5 cm  wide  notch  region  is  unperturbed  by  the  presence  of  the 
loading  pins. 

The  fracture  stress,  and  therefore  the  average  strain  energy  density, 
was  controlled  by  varying  the  depth  (aQ)  of  the  constant  radius  notch 
root.  When  a crack  erupts  from  the  notch  root,  it  is  subjected  to  higher 
strain  energies  than  would  be  expected  for  a sharp  crack  initiating  under 
a similar  loading  geometry.  The  situation  may  be  viewed  as  an  instan- 
taneous elevation  of  the  stress  field  intensity.  It  is  assumed  that  the 
apparent  stress  intensity  factor  for  this  initial  crack  is  reasonably 
described  by  the  stress  intensity  factor  characteristic  of  a sharp  crack 
having  length  equal  to  the  notch  depth,  but  the  situation  is  clearly 
dynamic  and  requires  calculation  of  an  appropriately  corrected  stress 
intensity  factor  as  given  by: 

dynamic  <a’V*o)  * ^static  (a’°)Hf(v)}  (2) 

where  Kgtat^c  (a*<0  is  the  static  stress  i-  isity  factor  calculated  for 
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f; 

FIGURE  4 


Modified  Single-Edge-Notch  Tension  Specimen  Used  to 
Investigate  Crack  Branching.  The  length-to-width  ratio 
(L/W)  in  the  gage  section  was  greater  than  4,  and  strain 
energy  density  was  varied  by  changing  the  notch  depth 
(a  ) . 
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a crack  length,  a,  in  a particular  specimen  geometry  loaded  by  stress, 
o.  The  factor  f(v)  is  the  dynamic  correction  to  the  stress  intensity 
factor  necessitated  by  the  effects  of  material  inertia.  According  to 
published  works  on  elastodynamic  cracks,  f(v)  is  a function  only  of 
crack  velocity. 

To  calculate  the  static  stress  intensity  factor  for  a sharp  crack 
that  has  catastrophically  initiated  from  a blunt  notch,  assume  that 
this  stress  intensity  factor  takes  the  following  form 

Vatu  <Vo)  3 ,KQ  (Vo)>  «> 

where  Kq  is  the  apparent  stress  intensity  factor  at  initial  fracture 
and  g is  a correction  factor  to  accommodate  extended  cracks,  a^  is  the 
instantaneous  crack  length,  aQ  is  the  starting  crack  length  (notch 
depth),  and  a is  the  fracture  stress. 

Among  the  solutions  for  crack  tip  stress  fields,  that  of  Tada, 
Paris,  and  Irwin  (10)  appears  to  be  the  most  precise  expression  for 
stress  intensity  factor  in  the  single-edge-notch  tension  specimen.  His 
result  is 

K - o/S-  . /(*.»  »/X  . °,.?52.  + 1-286X  + si,  X)3  (4) 

Q o cos  X 

where  X = ira  /2W 
o 

and  a is  applied  stress;  aQ  is  initial  crack  length  or,  for  present 

purposes,  notch  depth;  W is  specimen  width. 

This  expression  must  be  modified  to  represent  reasonably  the 
change  in  static  stress  intensity  factor  that  occurs  when  the  crack 
length  increases  from  a0  to  a^.  Because  the  specimen  geometry  was  de- 
signed with  particular  attention  to  eliminating  the  possibility  of 
stress  wave  reflections  from  the  boundaries  during  short,  rapid  crack 
runs,  it  is  reasonable  to  seek  a first  correction  to  the  solution  which 
describes  the  crack  as  if  it  had  statically  extended  in  an  infinite 
medium.  The  simplest  possibility  is  the  Griffith  crack,  for  which  the 
stress  intensity  factor  increases  as  the  square  root  of  the  crack 
length.  For  such  a flaw  the  stress  intensity  factors  for  cracks  of 
different  lengths  have  the  ratio 


Neglecting  dynamic  corrections  for  the  moment,  the  case  of  rapid  incre- 
mental extension  of  a crack  in  the  SEN  specimen  is  analogous  to  exten- 
sion of  a Griffith  flaw  in  an  infinite  specimen  for  purposes  of  calcu- 
lating static  stress  intensity  factors.  Therefore,  equation  (5)  is 
adopted  as  the  correction  factor  necessary  to  calculate  the  static 
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stress  intensity  factor  for  instantaneous  crack  length, 
expression  for  {equation  (4)}.  Specifically, 


ai» 


from  Tada's 


(6) 


The  velocity  correction  used  in  this  work  assumes  the  crack  is 
traveling  at  40%  of  the  Rayleigh  wave  velocity  (Cr) , which  is  assumed 
to  be  the  limiting  velocity.  This  is  based  on  several  considerations, 
the  first  of  which  is  the  experimental  work  of  Schardin  (11) , who  demon- 
strated with  high  speed  photography  that  the  speed  of  advance  of  the 
locus  of  crack  fronts  in  a process  of  fracture  by  multiple  branching 
was  on  the  order  of  95%  of  the  terminal  velocity  for  individual  cracks 
in  glass  specimens.  He  also  shows  that  increasing  the  applied  tension 
decreases  the  distance  to  initiate  a branch  point.  This  inverse  rela- 
tionship between  stress  and  crack  length  is  suggestive  of  the  stress 
intensity  factor  and  qualitatively  supports  the  notion  that 
Kg  ^ <j  (irag)  1/2  for  an  infinite  plate  may  be  a material  constant. 

The  observation  on  terminal  velocity  implies  that  the  crack  spends  very 
little  of  its  travel  time  undergoing  intense  acceleration  and  thus 
limiting  velocity  is  approached  in  a very  short  distance.  Further 
support  for  this  view  is  provided  by  Kerkhof  (12) , who  used  ultrasonic 
fractography  to  measure  crack  velocity  in  glass.  By  introducing  ultra- 
sonic waves  into  a tensile  fracture  specimen,  the  maximum  principal 
stresses  that  govern  the  direction  of  crack  propagation  are  made  to 
vary  periodically  in  space  and  time.  As  a consequence  the  fracture 
surface  topography  exhibits  waviness,  and  the  wavelength  is  determined 
by  the  ratio  of  the  crack  speed  to  the  frequency  of  the  ultrasonic 
modulator.  Very  precise  local  determinations  of  crack  speed  are  possi- 
ble using  simple  optical  microscopy  and  appropriate  lighting.  Kerkhof 
observed  that  the  crack  speed  was  approximately  constant  through  most 
of  the  mirror  zone  and  was,  in  fact,  the  terminal  velocity  at  the  point 
where  deep  surface  hackle  (microbranching)  began.  He  also  discusses 
why  surface  measurements  of  the  trace  velocity  of  non-uniformly  propa- 
gating cracks  can  lead  to  questionable  results.  A variety  of  common 
techniques  for  measuring  crack  speeds  has  been  described  elsewhere  (13). 
It  is  sufficient  to  note  that  no  presently  available  technique  for 
measuring  crack  speeds  in  metals  has  the  resolution  to  permit  conclusive 
statements  about  the  crack  speed  and  acceleration  in  the  pseudo-  mirror 
zone  of  brittle  metals  to  be  made.  It  is  here  that  attention  needs  to 
be  focused.  Consequently,  it  is  reasonable  to  assume,  by  analogy  to 
Kerkhof's  work,  that  the  best  estimate  of  crack  speed  at  incipient  bi- 
furcation is  locally  equal  to  the  terminal  velocity  in  brittle  steel 
(0.4  CR) . 

The  numerical  value  assigned  to  f(v)  in  equation  (2)  depends  upon 
which  crack  model  is  used.  Sih  and  Irwin  (14)  consider  the  problem  of 
evenly  spaced  radial  cracks  originating  at  a point  and  terminating  on 
a circular  locus  to  estimate  dynamic  crack-extension  force  (G) . The 
correction  factor  depends  upon  both  crack  speed  and  Poisson's  ratio  and 
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is  approximately  0.77  for  the  conditions  assumed  to  prevail  in  the 
present  work.  The  analysis  of  Freund  (15)  is  based  upon  a single 
dynamic  crack  and  implies  that  0.72  is  the  correction  factor.  These 
are  obviously  equivalent,  but  because  the  cracks  considered  here 
require  analysis  from  initiation  to  branching,  they  fall  more  naturally 
under  Freund's  analysis.  The  final  expression  adopted  for  calculating 
dynamic  stress  intensity  factors  in  the  single-edge-notch  tension 
specimen  may  now  be  written  as 

“dynamic  <ai’a0’a’v)  = °'72  /ai/ao  KQ  (a0'a)  (7) 

where  a^  is  the  instantaneous  length  of  the  dynamic  crack, 
a0  is  the  starting  notch  depth,  and  Kq  is  the  apparent 
stress  intensity  factor  given  by  equation  (4) . 

Using  this  formulation  of  the  dynamic  stress  intensity  factor,  it  is 
possible  to  estimate  the  strain  energy  release  rates  associated  with 
significant  features  detected  on  the  fracture  surface  during  post- 
mortem analysis  of  failed  specimens. 

Specimens  were  machined  to  the  geometry  shown  in  Figure  4 and 
heat-treated  to  the  conditions  listed  in  Table  1.  Crack  branching 
tests  were  performed  by  loading  the  specimens  to  failure  in  uniaxial 
tension  on  a universal  materials  testing  machine.  These  tests  were 
conducted  at  room  temperature  in  a displacement  control  mode  to  provide 
maximum  stiffness  in  the  load  train.  Quasistatic  displacements  were 
used  in  order  to  assume  "fixed-grip"  conditions  in  the  analysis.  This 
means  that  fracture  energies  come  only  from  the  high  strain  energy 
densities  stored  in  the  specimens  themselves  and  not  from  any  cross- 
head motion  in  the  testing  machine.  The  test  series  for  each  material 
condition  consisted  of  a range  of  notch  depths  selected  to  provide  a 
range  of  fracture  stresses. 

A typical  example  of  a full  branch  profile  is  shown  in  Figure  5. 
Such  samples  were  prepared  by  removing  unwanted  material  outside  the 
crack  branching  zone  with  an  abrasive  cut-off  wheel  and  reconstructing 
the  specimen  from  recovered  fragments.  The  assemblage  was  then  fixed 
in  Kold-Mount  and  a belt  sander  was  used  to  remove  half  the  specimen 
thickness.  Specimens  were  then  metallographically  polished  and  etched; 
the  view  in  the  photograph  is  a fracture  profile  along  the  vertical 
mid-plane  of  the  specimen.  Macroscopic  branch  cracks  were  found  which 
exhibited  all  the  behavior  patterns  alluded  to  by  Irwin  (16).  In 
general,  the  micro-branch  planes  were  oriented  at  angles  of  20-40°  to 
the  primary  crack  plane. 

Examination  of  the  fracture  surfaces  of  brittle  metals  revealed 
three  clear  features  that  are  analogous  to  features  found  on  glass 
fracture  surfaces  (Figure  6) . These  are  a pseudo-mirror  region  con- 
taining the  fracture  origin  centered  at  a distance  a^  from  the  specimen 
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FIGURE  5.  Crack  Branching  Profile  on  Mid-Plane 
of  FS-01  (B)  Steel  Specimen  (2.5  X). 
Note  the  skewing  and  subsequent  arrest 
of  the  main  crack. 


FIGURE  6 


Macroscopic  features  of  branching  crack  surface  in  a 
brittle  metallic  SEN  specimen.  Notch  depth  is  Sq? 
distance  to  fracture  origin  is  avj;  crack  length  at 
incipient  branching  (hackle)  is  aH;  and  crack  length 
at  full  bifurcation  is  ag. 
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edge,  a definable  houndary  to  that  region  beyond  which  deep  hackle 
occurs,  au,  and  a point  where  the  primary  fracture  fully  bifurcates,  aB. 
Hackle  ana  incipient  microbranching  are  considered  to  be  the  same. 
Because  the  boundary  separating  that  region  from  the  mirror  zone  is 
approximately  semi-circular,  ajj  was  defined  to  be  the  deepest  penetra- 
tion of  that  boundary  from  the  specimen  edge  at  which  the  fracture 
initiated.  The  detailed  results  of  the  crack  branching  tests  for  all 
materials  are  presented  in  Tables  6 through  12.  It  can  be  seen  that 
full  bifurcation  did  not  always  occur  even  in  the  most  brittle  condi- 
tions. Examination  of  the  data  for  materials  that  had  both  branched 
and  unbranched  fractures  indicates  that  insufficient  fracture  stress 
(or  strain  energy)  is  not  the  explanation.  It  was  observed  both  here 
and  in  the  work  of  Anthony  (4)  that  full  branching  is  the  culmination 
of  a build-up  of  microbranch  trials  until  one  of  these  microbranches 
can  develop  into  an  independent  crack.  Thus  there  appears  to  be  a 
statistical  quality  about  the  full  branch  crack  even  when  the  crack 
extension  forces  are  adequate. 

Because  hackle  always  occurred,  it  was  adopted  as  a more  convenient 
property  than  macrobranching.  The  dynamic  stress  intensity  factors  for 
hackle  (KH)  and  branching  (Kg)  are  independent  of  notch  depth  and  frac- 
ture stress  in  all  material  conditions.  This  fact  is  clearly  illus- 
trated by  Figures  7 through  10,  in  which  the  dotted  horizontal  lines 
are  the  calculated  mean  values  of  Kg  and  Kg.  Such  independence  lends 
further  credence  to  the  interpretation  of  Kg  and  Kg  as  material  prop- 
erties. 

SEM  Fractography  of  High  Energy  Fractures 

Numerous  fracture  surfaces  from  tensile  tests,  plane  strain  frac- 
ture toughness  tests,  crack  branching  tests,  and  explosive  fragmenta- 
tion tests  were  examined  by  scanning  electron  microscopy.  Typical 
results  are  illustrated  and  discussed  for  each  material.  All  four  con- 
ditions of  FS-01  steel  produced  fractures  that  were  essentially  the 
same  fractographically,  regardless  of  the  type  of  test  specimen  and, 
therefore,  the  magnitude  of  the  strain  rate  causing  specimen  failure. 
Figure  11  shows  the  basic  features  of  fracture  surfaces  of  FS-01  (C) 
steel.  This  material  is  a high  quality  tool  steel  with  a characteris- 
tically broad  dispersion  of  fine  carbides  that  serve  as  nuclei  for  void 
growth.  The  general  fracture  mode  seems  to  be  void -nucleated  quasi- 
cleavage. It  is  clear  from  the  figures  that  the  fracture  mechanism  is 
the  same  for  crack  branching  and  explosive  loading  test  specimens. 

Fractographs  of  HF-1  (A)  steel  are  shown  in  Figure  12.  Again  the 
fractures  are  identical  in  the  crack  branching  and  explosive  loading 
tests;  the  fracture  is  clearly  intergranular.  The  fracture  surfaces  of 
HF-1  (B)  steel  are  shown  in  Figure  13.  A high  velocity  fracture  sur- 
face near  a branch  point  from  a SEN  test  specimen  is  shown  in  Figure 
13A;  the  fracture  is  intergranular.  The  explosive  test  produced  the 
surface  shown  in  Figure  13B,  which  shows  additional  features.  The 
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FIGURE  7.  Dynamic  Stress  Intensity  Factors  for  Hackle  Versus 
Fracture  Stress  for  All  Conditions  of  FS-01  Steel. 


FIGURE  8.  Dynamic  Stress  Intensity  Factors  for  Crack  Branching 
Versus  Fracture  Stress  for  All  Conditions  of  FS-01 
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FIGURE  9.  Dynamic  Stress  Intensity  Factors  for  Hackle  Versus 
Fracture  Stress  for  All  Conditions  of  HF-1  Steel. 
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FIGURE  10.  Dynamic  Stress  Intensity  Factors  for  Crack  Branching 
Versus  Fracture  Stress  for  All  Conditions  of  HF-1 
Steel. 
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surface  has  been  somewhat  damaged  by  abrasion  and  oxidation  (white 
fringes).  There  are,  however,  undamaged  regions  that  show  the  same  in- 
tergranular fracture  as  the  branching  specimen.  The  most  striking 
feature  is  the  presence  of  shear  bands  which  are  the  thin  vertical  lines 
in  the  upper  right  hand  side  of  the  photograph.  What  is  significant  is 
that  the  microbranch  cuts  through  these  bands  which  are  in  approximate 
alignment  across  the  crack  opening.  This  indicates  that  the  bands  were 
there  first  and  the  brittle  radial  cracks  followed.  Figure  14  illus- 
trates the  intergranular  nature  of  the  fracture  in  HF-1  (C)  steel.  In 
spite  of  the  surface  damage  to  the  particle  recovered  from  an  explosive 
test,  the  primary  fracture  shows  a strong  intergranular  component.  The 
extent  of  grain  growth  in  the  HF-1  (A)  and  HF-1  (C)  heat  treatments  is 
also  evident  by  comparison  with  HF-1  (B) , which  had  little  or  no  grain 
growth  because  of  the  moderate  austenitizing  temperature. 

DISCUSSION  OF  EXPERIMENTAL  RESULTS  AND  COMPARISON  WITH  THEORY 


The  numerical  results  of  the  previous  sections  have  been  evaluated 
in  terms  of  mean  values  for  each  material  and  heat  treatment.  Standard 
deviations  have  also  been  computed.  Those  properties  relevant  to  the 
proposed  model  are:  plane  strain  fracture  toughness  (Kjc) , dynamic 
stress  intensity  factors  for  hackle  (Kg)  and  microbranching  (Kg) , and 
the  Mott  statistical  parameters  for  number  (N0)  and  size  (p) . A com- 
plete listing  is  provided  in  Table  13.  Calculations  indicated  that  a 
fourth  power  relationship  should  exist  between  the  Mott  number  and  size 
parameters  and  the  Kg/Kj  ratio.  Subsequent  experimental  results  indi- 
cated that  Kg/Kic  would  Be  more  appropriate  than  Kg/K^Q.  The  phenomen- 
ology of  crack  branching  itself  provides  a rationale  for  this  behavior. 
Prior  to  full  bifurcation,  the  microbranch  trials  produce  progressively 
longer  cracks,  the  lengths  of  which  are  greater  than  average  fragment 
dimensions.  It  is  for  this  reason  that  Kg  was  considered  to  be  not 
merely  more  convenient  but  perhaps  more  significant  than  Kg  in  the 
explosive  strain  rate  regime  in  which  brittle  metals  actually  shatter. 
However,  it  is  reasonable  that  the  crack  extension  force,  Gg,  for  a 
bifurcated  crack  (i.e.,  two  cracks)  should  be  about  twice  the  extension 
force,  Gg,  for  a single  dynamic  crack  about  to  attempt  branching.  The 
observed  ratio  Gg/Gg  was  found  to  be  1.7  (+0.5)  over  all  material  con- 
ditions. The  approximate  proportionality  of  Kg  to  Kg  further  justifies 
the  use  of  the  more  convenient  experimental  property,  namely  Kg. 

Figures  15  and  16  are  log-log  plots  of  the  Mott  size  parameter,  p,  as  a 
function  of  the  Kg/Kjc  ratio  and  the  Mott  number  parameter,  N0,  as  a 
function  of  the  Kg/Kj^  ratio.  The  best  straight  line  fitted  by  linear 
regression  through  the  HF-1  data  yielded  the  following  results: 


1.10  (Kg/Kjc)3’52  milligrams 

(8) 

2.82  x 105  (Kh/Kic)-3’ 30 

(9) 
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served  in  these  material  conditions. 


✓ 


FIGURE  15.  Mott  Size  Parameter,  y,  as  a Function  of  the 

VKic  Ra“°' 
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The  best  straight  line  fitted  by  linear  regression  through  the 
three  lowest-mass  points  in  the  FS-01  data  yielded  the  following 
results: 

pF  =0.14  (K^/K^)3'^  milligrams  (10) 

n/  = 3.16  x 106  (VK1C)"3*54  (11) 

The  exponents  in  equations  (8)  through  (11)  relating  the  Mott 
parameters  p and  Nc  to  the  Ky/fC^c  ratio  for  both  materials  are  equal 
within  experimental  uncertainties.  The  average  value  of  3.43  (+0.12) 
compares  favorably  with  the  model  prediction  of  4.  The  FS-01  (D)  data 
were  omitted  from  the  calculations  because  they  would  obviously  have 
distorted  the  results.  The  omission  is  justified  on  several  grounds. 
The  tests  barely  met  the  plane  strain  criterion,  i.e.,  fracture  speci- 
mens had  sizable  shear  lips;  moreover  the  particles  recovered  from  the 
explosive  tests  were  predominantly  bounded  by  shear  surfaces  and  there- 
fore not  suited  to  the  opening  mode  analysis  used  in  this  work.  That 
particular  datum  is  useful  in  that  it  sets  definite  bounds  on  the  model 
developed  here,  which  apparently  works  only  for  shatter  of  brittle 
cylinders  exhibiting  substantial  radial  components  of  fracture. 

A final  comment  is  in  order  regarding  the  nature  of  the  crack 
velocity  correction.  Both  experimental  limitations  and  analytical 
considerations  led  to  applying  the  same  correction  factor  to  all  static 
Ky  and  Kg  calculations.  The  finally  determined  functional  relation- 
ships are  simply  power  laws.  Omission  of  this  velocity  correction, 
i.e.,  using  the  static  values  of  stress  intensity  factors,  would  yield 
the  same  power  laws.  What  would  change  is  the  coefficients  rather  than 
the  exponents.  It  was  considered  important  to  have  these  coefficients 
reasonably  accurate  in  order  to  determine  whether  they  were  signifi- 
cantly influenced  by  metallurgical  factors  such  as  grain  size.  There 
are  indications  of  such  an  effect,  but  the  influence  of  grain  size  was 
not  systematically  investigated. 


SUMMARY  AND  CONCLUSIONS 


Dynamic  fracture  was  studied  in  two  hypereutectoid  steels,  FS-01 
and  HF-1,  that  easily  shatter  when  loaded  at  explosive  strain  rates. 
The  number  and  size  distribution  of  particles  generated  by  such  a 
process  is  well-characterized  by  the  Mott  cumulative  distribution 
function.  The  tendency  for  the  fragments  to  coarsen  with  increasing 
metal  ductility  is  a qualitative  observation  made  here  and  elsewhere. 
The  phenomenology  of  brittle  fragmentation  was  examined  and  a simple 
model  was  developed  based  on  the  concept  of  crack  branching  phenomen- 
ology. The  model  was  made  quantitative  by  establishing  relationships 
between  the  energies  associated  with  crack  branching  and  the  Mott 
parameters  that  characterize  the  particle  distribution.  This  model 
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predicts  that  the  average  mass  of  a particle  in  the  shattered  body  is 
proportional  to  the  fourth  power  of  the  Kg/Kj  ratio.  Kh/KIc  was  more 
convenient  experimentally  than  Kg/Kj^.  In  subsequent  experiments  the 
dependence  of  the  Mott  parameters  on  %/KIC  agreed  reasonably  well  with 
fourth  power  law,  so  long  as  there  was  a significant  amount  of  radial 
fracture  in  the  exploding  cylinder.  The  most  ductile  condition  of 
FS-01  steel  produced  results  that  were  far  out  of  line  with  predictions 
because  the  cylinders  fragmented  entirely  by  shear,  which  is  contrary 
to  the  assumption  of  brittle  fracture  in  the  model;  so  this  result  is 
not  surprising.  However,  it  is  noteworthy  that  static  tensile  elonga- 
tion correlated  monotonically  with  the  Mott  parameters  for  each 
material.  Since  these  high-strength  brittle  steels  are  not  particu- 
larly sensitive  to  strain  rate,  the  dependence  of  these  parameters  on 
the  ductility  may  be  through  a Kj  - ^Cod  “ ^ tyPe  functional, 
which  would  be  different  for  each  material.  This  is  believed  to  be 
the  first  time  that  dynamic  material  properties  were  used  to  construct 
a model  that  quantitatively  predicts  the  relative  size  distributions 
of  debris  in  blast- loaded  structures  of  different  materials. 


The  microscopic  fracture  modes  were  shown  to  be  essentially  inde- 
pendent of  strain  rate  throyghout, the  spectrum  from  quasistatic  to 


explosive  (i.e.,  10“^  to 
This  supports  the  fracture  mechanics  approach. 


throughout  the  spectrum 
:o  10°  sec*1) , at  least 


for  tensile  failures. 


The  energy  approach  used  in  this  study  may  be  more  general  than 
was  implied  by  the  restriction  to  brittle  opening  mode  cracks.  It 
should  be  possible,  although  much  more  difficult,  to  develop  similar 
models  for  Mode  II  fractures.  The  case  of  failure  by  full-wall  shear 
could  then  be  handled  and  the  mixed  mode  case  made  more  accurate. 


In  reality,  it  is  the  fracture  energy  that  dominates  the  break-up. 
The  proposed  crack  branching  modification  to  Mott's  fragmentation  model 
should  be  interpreted  symbolically,  rather  than  literally,  because 
parallel  cracks  in  close  proximity  to  each  other  can  propagate  simul- 
taneously only  if  they  are  both  traveling  at  terminal  velocity,  or  one 
will  be  arrested.  Normally  crack  branching  will  occur.  But  more 
importantly,  even  if  it  doesn't,  such  contiguous  cracks  will  experience 
stress  intensities  of  the  same  order  of  magnitude  as  for  microbranch- 
ing. Therefore  the  energies  will  be  the  same  as  for  the  multiple 
branching  events.  After  all,  two  very  closely  spaced  terminal  velocity 
cracks  cannot  yield  intrinsically  any  information  about  their  origins. 
Specifically,  they  cannot  indicate  whether  they  are  branches  of  a 
former  single  crack  or  whether  they  nucleated  independently.  This 
aspect  of  the  problem  suggests  a very  general  applicability  of  the 
method  of  analysis  to  all  problems  characterized  by  a multiplicity  of 
finely  spaced,  high-velocity  cracks. 

Besides  the  obvious  application  of  the  present  formalism  to 
military  fragmentation  devices,  it  may  be  of  value  in  studying  the 
inverse  problem  of  containment  of  accidental  explosions  without  the 
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formation  of  hazardous  high-energy  fragments.  Another  area  of  practi- 
cal technological  significance  to  which  this  approach  may  be  extended 
is  the  mechanics  of  rock  fragmentation.  By  characterizing  the  particle 
size  distributions  in  shattered  structures  of  common  rock  in  terms  of 
explosive  brisance  or  charge-to-mass  ratio,  it  may  be  possible  to 
specify  the  conditions  necessary  to  produce  a desired  crater  size  or 
to  produce  average  particle  sizes  which  may  be  conveniently  removed  by 
commercial  mining  equipment.  Along  the  same  lines,  the  present  formal- 
ism might  be  used  to  determine  the  most  efficient  impact  energies  for 
ore  pulverization  or  mineral  comminution.  There  are  apparently  many 
potential  applications  of  crack  branching  methodology  to  fragmentation 
of  brittle  materials  under  dynamic  loading  conditions. 
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